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The equation of state (EOS) of dense matter has been a long-sought goal of nuclear physics. 
Equations of state generate unique mass versus radius (M-R) relations for neutron stars, the ultra- 
dense remnants of stellar evolution. In this work, we determine the neutron star mass-radius relation 
and, based on recent observations of both transiently accreting and bursting sources, we show that 
the radius of a 1.4 solar mass neutron star lies between 10.4 and 12.9 km, independent of assumptions 
about the composition of the core. We show, for the first time, that these constraints remain valid 
upon removal from our sample of the most extreme transient sources or of the entire set of bursting 
sources; our constraints also apply even if deconfined quark matter exists in the neutron star core. 
Our results significantly constrain the dense matter EOS and are, furthermore, consistent with 
constraints from both heavy-ion collisions and theoretical studies of neutron matter. We predict a 
relatively weak dependence of the symmetry energy on the density and a value for the neutron skin 
thickness of lead which is less than 0.20 fm, results that are testable in forthcoming experiments. 

PACS numbers: 26.60.-c, 21.65. Cd, 26.60. Kp, 97.60.Jd 
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The masses of several neutron stars have been precisely 
measured using pulsar timing [l| ; simultaneous mass and 
radius measurements, however, are considerably less cer- 
tain. The leading candidates for such measurements are 
bursting neutron stars that show photospheric radius ex- 
pansion (PRE) 0] (for a review, see Lewin et al. 1993 H) 
and transiently accreting neutron stars in quiescence [J] . 

Observations have already begun to determine the uni- 
versal M-R relation for neutron stars and to place strong 
constraints on the EOS of dense matter. Previously de- 
rived constraints [H-Q have several limitations, includ- 
ing the use of fixed parametrizations for the EOS. Those 
works did not show their results to be independent of 
their parametrizations, including the possibility of de- 
confined quark matter. Neither did they address the full 
set of sources analyzed here. Ozel et al. [f| considered 
only PRE sources, but these may be subject to consid- 
erable systematic errors @, & EH ■ Steiner et al. Q con- 
sidered both types of sources but used a smaller data 
set. Here we use eight neutron stars: four which pro- 
duced PRE X-ray bursts (4U 1608-522 [TO], KS 1731- 
260 [T3, EXO 1745-248 [lj, and 4U 1820-30 Q) and 
four in quiescent low-mass X-ray binaries (qLMXBs) in 
the globular clusters M13 fll , u Cen [HI , 47 Tuc and 
NGC 6397 0. Refs. d|l| took advantage of modern 



predictions for pure neutron matter near the saturation 
density to constrain the M-R relation, but ignored sys- 
tematic uncertainties associated with the observations. 
In this paper, we demonstrate that the inferred M-R 
constraints are insensitive to the removal of either all the 
PRE burst sources or the most extreme qLMXB sources. 

At the lowest energy densities (< 15 MeV/fm 3 or 
< 3 x 10 13 g/cm 3 ) the pressure-density relation is well- 



understood pi]. Between 15 and 200-300 MeV/fm 3 , the 
EOS is well-described by four parameters, the incom- 
prcssibility, the skewness, the magnitude of the symme- 
try energy (S v ), and the parameter describing the den- 
sity derivative of the symmetry energy (L), all evalu- 
ated at the nuclear saturation density (approximately 
150 MeV fm -3 ). These parameters are constrained to 
varying degrees by experimental data [2(| , including nu- 
clear masses [2lj], neutron skin thicknesses I22I . giant 
dipolc resonances and dipole polarizabilities [23l425l | , and 
heavy- ion collisions [26jj. High-density matter is con- 
strained by (i) causality (the speed of sound must not 
exceed the speed of light), (ii) hydrodynamical stabil- 
ity, and (iii) having a sufficient maximum mass (it must 
be greater than the largest well-determined neutron star 
mass, 1.97 ± 0.04 M Q for PSR J1614-2230 [13). In addi- 
tion, we impose the constraint that implied neutron star 
masses cannot be less than what is achievable in super- 
nova, about 0.8 Mq. 

The low-density part of all the EOS parametriza- 
tions (except for strange quark stars) are described as 
in Stciner et al. 0- All of the parameters in the low- 
density EOS, as well as the parameters in the high- 
density EOS models, are described with uniform prior 
distributions. Our fiducial EOS model (A) parametrizes 
the high-density EOS as a set of two piecewise continuous 
power laws defining pressure P = e L+1 ' n as a function 
of energy density e. Model A has four high-density pa- 
rameters: the transition energy density between the low- 
density EOS and the first polytrope, the transition en- 
ergy density between the first and second polytrope, and 
the two poly tropic indices, n\ and n 2 . We also employ an 
EOS model (B) which is similar to model A, except that 
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the exponents in the two polytropes are parametrized 
with uniform priors in I\ instead of rij. A third EOS 
model (C), parametrizes the EOS at high densities with 
a uniform prior in the pressure at four fixed energy densi- 
ties, 400, 600, 1000, and 1400 McV/fm 3 . The low-density 
EOS is used up to energy densities of 200 MeV/fm 3 and 
the EOS is assumed to be linear between that point and 
400 MeV/fm 3 . The linear relation between energy densi- 
ties of 1000 and 1400 MeV/fm 3 is extrapolated to higher 
energy densities when necessary. All parameters are cho- 
sen with large enough ranges to ensure the results do not 
change significantly when the range is increased. 

A fourth EOS model (D) assumes that matter is a 
polytropc at intermediate densities and quark matter at 
high densities. The pressure of quark matter is described 
by [H| 
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where fi is the quark chemical potential. The quantity B 
is the bag constant and simulates confinement. The pa- 
rameter ci4 describes corrections to the leading coefficient 
from a non- interacting Fermi gas (for which 04 = 1). Cor- 
rections from pcrturbativc quantum chromodynamics at 
high density suggest 0.6 < 04 < 1. The parameter ai ap- 
proximately describes corrections from the finite strange 
quark mass m s and the quark supcrfluid gap A and is 
given by a 2 = m 2 s — 4A 2 (see also Ref. Q). In order to 
include the effects of a possible mixed phase in a model- 
independent way, a polytrope is added in between the 
low-density EOS and quark matter. In this model, there 
are five parameters in total: the transition energy density 
between low densities and the polytrope, the polytropic 
index, the transition energy density between the poly- 
trope and quark matter (used to fix the value of B), 02, 
and 04. To describe bare strange quark stars, Model E 
applies Equation 1 at all densities, with neither a low- 
density EOS nor an intermediate polytrope. Model E 
thus has three parameters, 02, 04, and B. 

Our baseline data set includes all eight astrophysical 
sources, interprets the PRE burst sources with an ex- 
tended photosphere, and assumes the same distribution 
of color correction factors and distances as in Steiner, et 
al. 0. We consider several modifications to the baseline 
scenario. Suleimanov et al. [Io[ have suggested that the 
X-ray spectra for PRE sources are affected by accretion 
and the eclipse of the neutron star by the disk. This 
affects the normalization at late times and we take this 
into account by increasing the color correction factor fc, 
taking 1.45 < fc < 1.8 (modification I). Alternatively, 
Boutloukos et al. have suggested that the color cor- 
rection factors are, instead, smaller as a result of mag- 
netic confinement of the X-ray burst; this is considered 
in modification II in which we assume 1 < fc < 1.35. 
Some previous works assumed the photosphere of PRE 
bursts is coincident with the neutron star surface [dl , and 
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11.23 


11.53 


12.17 


12.45 
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10.63 


10.88 


11.45 


11.83 
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11.44 


11.69 


12.27 


12.54 


Variations in the data interpretation 


A I 


11.82 


12.07 


12.62 


12.89 


A II 


10.42 


10.58 


11.09 


11.61 


A III 


10.74 


10.93 


11.46 


11.72 


A IV 


10.87 


11.19 


11.81 


12.13 


A V 


10.94 


11.25 


11.88 


12.22 


A VI 


11.23 


11.56 


12.23 


12.49 


Global limits 


10.42 


10.58 


12.62 


12.89 


More speculative scenarios 


C II 


9.17 


9.34 


9.78 


10.07 


A VII 


12.14 


12.29 


12.63 


12.81 


E 


10.19 


10.64 


11.57 


12.01 


A VIII 


12.35 


12.83 


13.61 


13.92 



TABLE I: Limits for the radius of a 1.4 solar mass neutron 
star for all of the models considered in this work. Model A 
and the assumption 1.33 < fc < 1-47 for the PRE sources is 
assumed unless specified otherwise. 



we also consider this scenario (III). Finally, we test the 
sensitivity of our results to the removal of any one source 
or class of sources by removing X7 (IV) or M13 (V), and 
by removing all PRE sources (VI). To summarize, the 
standard models we examine include the baseline case 
(A), three variations of the EOS (B-D), and six other 
modifications of the baseline model varying the included 
data or its interpretation (A I-A VI) . In addition, we also 
examine several more speculative scenarios which are de- 
scribed below. 

We use the Bayesian method of Ref. Q , using marginal 
estimation to determine the posterior probability densi- 
ties of quantities of interest. The marginal estimation in- 
tegrals are performed using Markov Chain Monte Carlo. 

In our baseline analysis, we find strong constraints on 
the M-R curve and on the dense matter EOS: the radius 
of a 1.4 M Q neutron star is between 11.2 and 12.3 km 
(95% confidence). The permissible radius range encom- 
passing all variations of the EOS and interpretations of 
the astrophysical data, but not including the more spec- 
ulative scenarios, is 10.4-12.9 km (95% confidence), only 
moderately larger than the baseline result. The 68% and 
95% confidence ranges are displayed in the upper and 
middle portions of Table |U We determine the M-R rela- 
tion for a range of neutron star masses (see Figure 1). We 
also determine the 68% and 95% confidence intervals of 
the EOS of dense matter (Figure 2). The estimated un- 
certainty of the pressure is approximately 30-50% at all 
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FIG. 1: A comparison of the predicted M-R relation with the 
observations. The shaded regions outline the 68% and 95% 
confidences for the M-R relation; these include variations in 
the EOS model and the modifications to the data set (see 
Table [TJ> but not the more speculative scenarios. The lines 
give the 95% confidence regions for the eight neutron stars in 
our data set. 



FIG. 2: The predicted pressure as a function of baryon den- 
sity of neutron-star matter as obtained from astrophysical 
observations. The region labeled "NS 68%" gives the 68% 
confidence limits and the region labeled "NS 95 %" gives the 
95% confidence limits. Results for neutron-star matter from 
effective field theory [jlj (see inset), from quantum Monte 
Carlo 13011 . and from constraints inferred from heavy-ion col- 
lisions [3l| are also shown for comparison. 



densities achievable in neutron star interiors. In addition, 
the posterior distributions of the central energy density 
of the maximum mass star imply that the highest central 
density is - 1200 MeV/fm 3 

Producing significantly different neutron star radii re- 
quires extreme assumptions regarding the EOS and the 
data. We now consider more speculative scenarios, which 
are presented in the bottom portion of Table U (see also 
Figure 3). To achieve significantly smaller radii, we must 
assume both that the color correction factor is anoma- 
lously small (< 1.3) for all of the PRE sources and that 
the EOS has strong phase transitions (model C). In this 
case, we get radii as small as 9 km. Increasing the max- 
imum mass constraint, as would be the case if the es- 
timated most-likely mass of the pulsar B1957+20 is 2.4 
M Q [32[, slightly increases radii (modification VII). We 
obtain even larger radii if we add the long PRE burst 
source 4U 1724-307 |10f and further assume, as suggested 
in Suleimanov et al. [lfj (modification VIII), that the 
short PRE bursts and the qLMXBs M13 and u Ccn not 
be considered because of modifications to their spectra 
due to accretion. This scenario cannot yet explain, how- 
ever, why short PRE burst cooling tails are observed to 
have constant normalizations. 

While we are able to significantly constrain the P-e re- 
lation, determination of the composition of neutron star 
cores is not yet possible. To this end, we consider EOS 



model E, which describes the entire star by the high- 
density quark matter EOS used in model D, i.e a self- 
bound strange quark star. In the mass range 1.4-2 solar 
masses, the radii are not significantly different from our 
baseline model so that there is no strong preference for 
cither strange quark or hadronic stars; however, model E 
predicts radii significantly less than 10 km for low masses 
(< 1.2 M ). 

Our neglect of rotation is unlikely to affect our conclu- 
sions. Rotation increases the radius at the equator and 
decreases the radius at the poles, and this could be rele- 
vant for the interpretation of some PRE X-ray bursts: the 
rotation rate of 4U 1608-522 is 619 Hz, more than half 
of the rate for which the equatorial radius is increased by 
about 50%. However, this is likely to produce a system- 
atic uncertainty smaller than that due to variations in 
fct which we have already taken into account. The rota- 
tion rates for the qLMXBs in our sample are unknown. 
Assuming they are similar to other qLMXBs, however, 
means that the effect of rotation is smaller than that of 
their distance uncertainties. 

The relationship between pressure and energy density 
(Figure 2) that we determine from our baseline analysis 
from observations is consistent with effective field the- 
ory [lj| and quantum Monte Carlo [U, [3(| calculations 
of low-density neutron matter. Note that these neutron 
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FIG. 3: Predicted M-R relations for different EOS models 
and data interpretations. Proceeding from back to front, 
the red contours and probability distributions are for strange 
quark stars (EOS model E with no modifications to the data) . 
Next are green contours which correspond to the baseline 
model (EOS model A with no modifications to the data set), 
and the magenta results are those assuming a larger maximum 
mass to accomodate a mass of 2.4 solar masses for B1957+20. 
Finally, the black lines are the 10 Skyrme models from Stone 
et al. which are inconsistent with the data because their radii 
are too large (they do not match the observations at masses 
low enough to accomodate the low-mass objects like M13). 



matter results are incompatible with Sulcimanov's inter- 
pretation of 4U 1724-307 and suggested exclusion 
of short PRE bursts and qLMXLBs M13 and w Ccn. 
Our results are also consistent with the high-density con- 
straints from heavy- ion collisions [3l| . 

Our results imply that over one third of the modern 
Skyrme models studied in Stone et al. (33j are inconsis- 
tent with observations. Covariant ficld-thcorctical mod- 
els that have symmetry energies which increase nearly 
linearly with density, such as the model NL3 34[, are 
also inconsistent with our results, although they may still 
adequately describe isospin-symmetric matter in nuclei. 

Our models do not place effective constraints on the 
symmetry parameter S v , but do place significant con- 
straints on the symmetry energy parameter L; these are 
summarized in Figure 4. The probability distribution for 
each model is renormalized to fix the maximum probabil- 
ity at unity and are then shifted upwards by an arbitrary 
amount. The range which encloses all of the models and 
modifications to the data is 43.3 to 66.5 MeV to 68% 
confidence and 41.1 to 83.4 MeV to 95% confidence. The 
allowed values of L arc substantially larger for Model 
C, which allows strong phase transitions at high densi- 



FIG. 4: The limits on the density derivative of the symmetry 
energy, L. The single-hatched (red) regions show the 95% 
confidence limits and the double-hatched (green) regions show 
the 68% confidence limits. 



ties, because the parametrization decouples the low- and 
high-density behaviors. 

Our preferred range for L is similar to that obtained 
from other astrophysical studies _ 
imcntal studies, e.g., Refs. [24, 




and exper- 
Our results sug- 
gest that the neutron skin thickness of 208 Pb [3(1, [37J 
is less than about 0.20 fm. This is compatible with ex- 
periment 38[ and also with measurements of the dipolc 
polarizability of 208 Pb 0. 

While we have endeavored to take into account some 
systematic uncertainties in our analysis, we cannot rule 
out corrections due to the small number of sources and 
to possible drastic modifications of the current under- 
standing of low-mass X-ray binaries. Nevertheless, it is 
encouraging that these astrophysical considerations agree 
not only with nuclear physics experiments but also with 
theoretical studies of neutron matter at low densities and 
heavy-ion experiments at higher densities. 
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